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Abstract 

Let d = (di, d,2, ■ ■ ■ , d n ) be a vector of nonnegative integers. We study the number 
of symmetric 0-1 matrices whose row sum vector equals d. While previous work has 
focussed on the case of zero diagonal, we allow diagonal entries to equal 1. Specifically, 
for D £ {1,2} we define the set Qo{d) of all n x n symmetric 0-1 matrices with row 
sums given by d, where each diagonal entry is multiplied by D when forming the row 
sum. We obtain asymptotically precise formulae for |£/D(tZ)| m the sparse range (where, 
roughly, the maximum row sum is o(n 1 / 2 )), and in the dense range (where, roughly, the 
average row sum is proportional to n and the row sums do not vary greatly). The case 
D = 1 corresponds to enumeration by the usual row sum of matrices. The case D = 2 
corresponds to enumeration by degree sequence of undirected graphs with loops but 
no repeated edges, due to the convention that a loop contributes 2 to the degree of its 
incident vertex. We also analyse the distribution of the trace of a random element of 
Qi)(d), and prove that it is well approximated by a binomial distribution in the dense 
range, and by a Poisson binomial distribution in the sparse range. 



1 Introduction 

Let d = (c?i, g? 2 , . . . , d n ) be a vector of nonnegative integers. Define G(d) to be the number 
ofnxn symmetric matrices over {0, 1} with zero diagonal, such that row j sums to dj, for 
j = l,...,n. 
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The quantity G(d) has been well studied, as cited below. In this paper we consider the 
case where the diagonal need not be zero. For D 6 {1,2} define Qo(d) to be the set of n x n 
symmetric matrices A = (a^) over {0, 1} such that 

D a jj + a jk = dj for j = 1, . . . , n. 

l<k<n, k^j 

We wish to find an asymptotic formula for 

G D (d) = \g D (d)\. 

The case of D — 1 corresponds to enumeration by row sum of symmetric 0-1 matrices. 
If we interpret A as the adjacency matrix of a simple undirected graph with loops, then 
the case of D = 2 corresponds to enumeration by degree sequence of simple undirected 
graphs with loops. Such graphs arise in various applications including the study of graph 
homomorphisms [9] and sign patterns jl]. 

Throughout the paper we will refer to a nonzero entry on the diagonal of a 0-1 matrix as 
a loop. For £ — 0, 1, . . . , n, let Go(d, £) be the set of matrices in Gn{d) with exactly £ loops 
(that is, with trace £), and let Gf)(d,£) = \Q D (d,£)\. Clearly we have Gi>(d, 0) = G(d) and 
Go{d) = X^=o Go(d, £). We also note here that G\(d,£) = unless ^2™ =1 dj has the same 
parity as £, and G2^d, £) = unless 5^? =1 dj is even. 

When dj = d for j = 1, . . . , n, we write Go(d) = GD{n,d) and refer to this as the regular 
case. 

We will use the following parameters frequently: 

n g 

A = -, d max = max dj, 

n — 1 j 

n n 

R=J2 ( d J - d Yi S r = J2 [ d A r ^ = 2 ' 3 )' 

3=1 3=1 

where [a] r = a(a — 1) ■ • ■ (a — r + 1) denotes the falling factorial. 

Throughout the paper, the asymptotic notation 0(f{m)) refers to the passage of the 
variable m to infinity. (Usually m = n or m = S.) In the dense setting we also use a 
modified notation 0(/(n)), which is to be taken as a shorthand for 0(/(n)n ce ) with c a 
numerical constant (perhaps a different constant for each occurrence). We write Q(g(n)) to 
indicate any function which is greater than Cg(n) for some constant C > and sufficiently 
large n. 
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It appears that there is very little prior research on Gn{d). The most general result, by 
Bender and Canfield, dates from 1978. 

Theorem 1.1 Q3J). Suppose that 1 < c? max = 0(1). Then 

Gi(d) = 7=2 (T(W'TA^-i - f - # + 

uniformly as S — > oo. 

Note that G\(n, 1) is the number of involutions on n letters (and also the number of 
Young tableaux with n cells, see [20J A000085]). The asymptotic expansion of G\(n, 1) was 
previously known, see [5l[T9] . We found no prior asymptotic work on G%(d) at all. 

In the case of D = 1, a graph with n vertices and I loops can be mapped to a graph with 
n + 1 vertices and no loops, by introducing a new vertex and replacing each loop by an edge 
to this vertex. This mapping is bijective and hence 

G 1 ((d 1 ,...,d n ),£) = G((d 1 ,...,d n ,£)). 

However, this doesn't seem to be of much use in asymptotic enumeration, since the important 
values of £ place the degree sequence (c?i, . . . , d n , i) out of range of existing explicit estimates. 

Our approach to estimating Gjj{d) will be to sum over all possible diagonals using the 
existing estimates for G(d). The main estimates we will use are the following two theorems. 
The history of previous results on G(d) is summarized in [14] and |16j . 

McKay and Wormald [171 Theorem 5.2] proved the following asymptotic formula for G(d) 
in the sparse regime. 

Theorem 1.2 ([17]). If 1 < d max = o(S^ 3 ) then 

S\ 

° {d) = (5/2)12^^=1^-! 6X 
uniformly as S —> oo, with S even. 

In the case of dense matrices, the following result was due to McKay and Wormald [16] 
except that we will use an improved error term from a generalization by McKay [15J. A 
less explicit formula allowing a wider variation of the degrees was proved by Barvinok and 
Hartigan [2]. 



C C2 C C4 C2 / J.i 
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Theorem 1.3 Q15J). Let a, b > be constants such that a + b < |. Then there is a constant 
£o = £o(a,b) > such that the following holds. Suppose that dj — d is uniformly O(n l l 2+£0 ) 
for j — 1, . . . , n and that 

n 

min{d, n — d — 1} > — 

3a log n 

for sufficiently large n. Then provided S is even we have 

<m = V2 - A)-)® exp( I - w{1 R ' x)2ni + 0(n->)) f[ (1.1) 

This formula also matches the sparse case under slightly more restricted conditions than 
Theorem II .21 and is conjectured to hold in the intermediate domain as well (see [T8| Theorem 
2.5] and the conjecture stated immediately thereafter). 

Note that Theorem 11.31 remains true if e (a,b) is decreased (but is still positive), since 
the conditions of the theorem become stronger. 

We now state our main enumeration theorems, starting with the dense regime. 

Theorem 1.4. Let a, b > be constants such that a + b < |. Then there is a constant 
e = e(a, b) > such that the following holds. Suppose that dj — d is uniformly 0(n l l 2+£ ) for 
j = 1, . . . , n and that 

n 

mm\d,n-d] > — (1.2) 

1 ' 1 ~ 3alogn v ' 

for sufficiently large n. For D G {1, 2}, define 

d 



n + D-1 



and let 

Qi(d r 



1 {t-df {1-dfR R 2 



4 Ad(n-d) 2d 2 {n-d) 2 Ad 2 {n-d) 2 ' 



Q 2 {d,£) = -- — — -=■ - —j- —^-7 + 



4 /i 2 (l - ^-i)n 2 4/x|(l - /i 2 ) 2 n 4 /i 2 (l - ^n 2 
(1 -2fi 2 )(e-fi 2 n)R 2(£-fx 2 n) 2 R 
+ ^(l-/x 2 )2n3 ^(l-^ 2 )% 4 ' 

When I has the same parity as S we have 



G x {dJ) = V2 (tf (1 - ^f' 2 (f)^ (1 " » 



(n-£)/2 



exp(Q 1 (d,£) + 0(n- b ))H 

3=1 
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while for £ = 0, . . . , n and even S we have 

c 2 (d, i) = v2 (tf (i - n 2 f-^) Q a (i - ^ 



e^{Q 2 (dA) + 0(n- b ))\["' 1 



j'=i 



(1 j 



Defining 



we have 



d x l 2 n - dn 



dV2 + (n-d)V2' ^ 2 ^ n+1 



= " (i - ah) 1 -*)" /2 (ni /2 +(i -/iO i/2 r 

n 

xexp(Q 1 (d,£ 1 ) + 0(n- b )) J] 
and, /or even 5, 

G 2 (d) = v/2 0<(1 - Ai2 ) 1 -^)(" 1 )exp(g 2 (d,£ 2 ) + 0(n- & )) f[ (^j. □ 

In Theorem 11.61 we will prove that £d is close to the expected number of loops in a 
randomly chosen element of Qo{d). For the reader's convenience, we note that 

Q 2 (d,£ 2 ) = -Ul R U - " 



A\ /z 2 (l - /i 2 )n 2 / V /i 2 (l -/i 2 )n 2 j 
Unfortunately, the expression for Qi(d,£i) does not simplify much. In the case of regular 
graphs we have R = 0, so the formulae for Q£>(d,£) simplify greatly and in particular 

— 71 

Qi(d,h)- 



2n + 4s/d(n - d) ' 
Our main result for the sparse case is the following. 
Theorem 1.5. Suppose that 1 < d max = o^S 1 ^ 3 ). Then 

i S/ VA ..V 1 f n: 1 S* Si 



4 S AS 2 



1 Q Q Q2 c2 c c4 o2 / J3 

( Do Dr. DnO'i Dr. On „ / O m 



24S 1 ' 2 S 3 / 2 3S 3 / 2 2S 5 / 2 2S 4 AS 5 6S 3 V 5 
uniformly as S — > 00 , and 

C2 C2c C4 C2 /J3 



G 2 (d) = v^(f) 5/2 (nd J !)" 1 exp(| 

uniformly as S —> 00 S 1 even. 



4S 2 2S 4 AS 5 QS 3 \ S 
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If S is even then we may replace the factor ^2 (S/e) s ' 2 by S\/((S/2)\ 2 s ' 2 ) . In the regular 
case the formulae simplify as follows. 

Corollary 1.1. Suppose that 1 < d — o(n 1 / 2 ). Then 

1 fnd\ nd/2 
G l {n,d) = -j={dr n 

(2-2d-d 2 24(n- l)d + 20d 2 + 11 d 3 I ' d 2 

x exp 1 — = — — + O — 

V 4 2AV^d 12n V™ 

uniformly as n — > oo, and 



( nd\ 1 f 
G 2 (n,d) = V2(d\)- n (—J expf 



(d— 3) rf 3 + C Y^_ 



4 12n \ n 



uniformly as n — >■ oo wift no? even. 



Again, if n<i is even then the factor \/2 (nd/e) nd ^ 2 may be replaced by (nd)!/((nd/2)! 2 nd / 2 ). 

Theorems 11.41 and 11.51 are proved in Section [2] and [31 respectively. Along the way we 
prove some technical results (Lemmas 12.11 13.21 13. 3ft which may be of independent interest. 
But first, in Section 11.11 we state a theorem on the distribution of the trace of a random 
element of Go{d), and discuss some interesting features of this distribution. Theorem 11.61 
is proved in Section HI Finally in Section [5] we state a conjecture regarding the number of 
regular graphs with loops, for all possible degrees. 



1.1 The distribution of the trace 

The calculations we will give in the process of proving Theorems ll.4l and fT31 will provide some 
information on the distribution of the trace of a random element of Qn{d)- We summarize 
that information here. 

For p = (pi, . . . ,p n ) G [0, l] n , let Xi, . . . ,X n be independent random variables with 
Prob(Xj = 0) = 1 — pj and Prob(Xj = 1) = pj for each j. The Poisson binomial distribution 
PB(p) is the distribution of Y^j=i-^-j- Define 

n 

PB(p, t) = Prob ( X i = £ ) ■ 

3=1 

The special case p = (p,p, . . . ,p) gives the familiar binomial distribution, 
PB((p, ...,p),£) = Bin(n, p, £) = (^V (1 " vY' 1 - 
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Theorem 1.6. Let Yd = Y^{d) be the random variable given by the trace of an element of 
Qo{d) chosen uniformly at random. 



(i) If the conditions of Theorem \l.J\ hold then, for I — 0, . . . , n, 

Prob(Yi = £) = (2 + 0(n- b ))Bin(n,l 1 /n,£) + 0(e- nn(1) ), 
E(Y 1 )=l l (l + 0(n- b )), 
Var(yi) = 4(1 - £i/n)(l + 0(rT 6 )), 

Prob(F 2 = £) = (1 + 0{n- b )) Bm{n,I 2 /n,£) + 0{e- nSm ), 
E(Y 2 )=£ 2 (l + 0(n- b )), 
Var(F 2 ) = l 2 {\ - I 2 /n){\ + 0(n- b )), 

where £ must have the same parity as S in the D = 1 case and S must be even in the 
D = 2 case. 



(ii) Define p' = (p[, . . . ,p' n ) and p" = (p'(, . . where for j = 1, . . .,n, 

Pj or- "f 

Pj 



dj d 3 {2d 3 -l) d] d j (d 3 -2)S 2 djSj 
y/S 2S S 3 / 2 S 5 / 2 2SV 2 ' 

j, d {d 3 -l) 



S 

If the conditions of Theorem \1.5\ hold then, for £ = 0, . . . , n, 

'dl 



Prob(F! = £) = 


v 2 + °( 


Eft) = Vs - ^ 


o c 

Var(Yi) = VS - -| 


Prob(F 2 = £) = 




E(y 2 ) = 




Var(F 2 ) = 





s 



+ S- 1 ^ PB(p',l) + 0{e 



1 fd 3 \ 
Of max • 



2 \ >5 1/2 / 



2 1 + 



d 3 ^ W S- 



max 



'2 



s jj s 

dl„A\ s 2 



S JJ S' 

where £ must have the same parity as S in the D = 1 case and S must be even in the 
D = 2 case. 
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The parameter jj, D can be thought of as measuring the density of entries equal to 1, 
while Yo/n is the density of loops in a randomly chosen element of Qo{d). In the dense 
range of Theorem 11.61 we see that Y 2 /n is concentrated near the same value [i 2 , while Y\ /n 
is concentrated near 

ii _ ygi 

Figure [T7T1 illustrates this curious difference between D = 1 and D = 2. 




o¥ — . , ' 1 . 1 . , ' 1 

0.2 0.4 0.6 0.8 1 

Figure 1: The expected density of the diagonal as a function of the overall density \ib- 

When D — 1, Theorem 11.61 tells us that the most significant term in E(Yd) depends only 
on 5* and not on d, within the range of d values allowed by the theorem. To explore this 
further, let A n (S) be the set of all n x n symmetric 0-1 matrices with exactly S entries equal 
to 1. The number of matrices in A n {S) with exactly I loops is 




when S and I have the same parity, and otherwise. For 1 < S < n 2 — 1, it can be proved 
that the maximum value of this function occurs either at i\ rounded up to an integer of the 
same parity as S or i\ rounded down to such an integer. On the basis of experiments, we 
conjecture that the mean number of loops in A n (S) always lies in (Ji — + |). 

Also note that I\ ~ \fS for 5* = o(n 2 ), matching the leading term of E(Yi) in the sparse 
case. 



S 



When D = 2 we consider instead the set B n (S) of all graphs with loops allowed, with 
n vertices and S/2 edges (loops counting twice). Matrices which correspond to graphs in 
B n (S) can be formed by choosing S/2 entries on or below the main diagonal, setting these 
equal to 1, then adding this matrix to its transpose. (Nonzero entries on the diagonal all 
equal 2, which is their contribution to the row sum.) The number of graphs in B n (S) with 
exactly I loops is 

Af (I) 
I J \S/2 - 

Up to scaling, this is the hypergeometric distribution with parameters (("^ )> n, S/2) and 
mean S/ (n + 1) = \xin. 

The binomial distributions in part (i) of the theorem are asymptotically normal, as is well 
known. The Poisson binomial distributions in part (ii) of the theorem are asymptotically 
normal for Y\ (see [8J), and asymptotically Poisson for Y 2 , by Le Cam's Theorem [TT] (see 
also [1, Equation 1.1]). 

2 The dense case 

In this section we prove Theorem 11.41 
2.1 A technical lemma 

We will require a technical lemma which might be of some independent interest. If f3 = 
. . . , (3 n ) is a vector of real numbers and I — 0, . . . , n, define 

e 

l<j\<—<ji<n s=l 

Lemma 2.1. Define (3 = \Y^j=xPj an ^ su PP ose that f3j — (3 = 0(n -1 / 2 ) uniformly ft 
j = 1, . . . , n. Then, for sufficiently small e > 0, we have 



or 



r- 

uniformly for t = 0, . . . , n. 



Proof. The factor e e/3 can be removed by replacing each by (3j — f3, so it suffices to prove 
the lemma for Y^j=i Pj = P = 0- 
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We divide the proof into three parts, depending on I. Let B = maxj\(3j\. Choose a 
constant c > such that Bn}l 2 ~ ce = o(l). 



First assume that n 1 ^ 2 C£ < £ < n — n l l 2 C£ . Since Ug(f3) is the coefficient of y l i 



m 



rL=i(l + e v)i we can estimate it using the saddle point method. We choose the contour 



to be a circle of radius r centered at the origin, where 



n-l 

For j — 1, . . . , n let 

~ i + e Pj r ' 

Changing variable according to y = re %9 and applying Cauchy's theorem, we obtain 

Um = P((3) r F(6) d9, 

J — n 

where 

p(a) n- =1 (i+^o n? =1 (i+^-i)) 

The coefficient ipj satisfies 



= -(l + 0(n- 1 / 2 )). 

We now divide the domain of integration into the two subdomains \9\ < 9 and \9\ > 9 , 
where 

Expanding F(6) for \9\ < 9 , we find using (12. ip that 

(n n 

n 

\ 2n n / 



where the 0(1) term is independent of 9. Since the interval \9\ < #o is symmetric about 0, 
we can instead integrate 

1 + = exp (-^^ # 2 + 0{n-^ . 
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Furthermore, 

j\ w (^JJ^_A 0i + din- 1 ^ d9 = n ~ n ^ 
(and similarly for the lower tail) and hence 



/>>^^-p(o(^,). 

For the complementary subdomain \9\ > 9 , note that 

|l + ^.( e * _1)| = v /l-2^(l-^)(l-cos^), 

which is a decreasing function for 9 G (O ,ir). Therefore, < |-F(#o)| f° r \@\ > $o- Since 

1 — cosy > 2y 2 /n 2 when —ix < y < it, we have 



Wo) | = II ^-^(l-^Oll-cos^o) 

( 21og 2 n ~ _ 1/2 A 
expl — hO(n 1 ) 1 



< 

-f2(logn) 



= n 



Hence 



F(<0 dB = + £° F(9) d9 = yQ^J ex P (0(n- 1 / 2) ). 



Finally, we calculate that 



n 



p(3) = ? n 1 + — 

Therefore 

^(0) = ^! exp f V /3 2 + O^ 1 / 2 )^) , 

which equals the expression in the lemma, by Stirling's formula. 

We next consider the case that < £ < n 1/2 ~ C£ . Expand U e ((3) = J2 s >o T s/ s ^ where 

l<ji<—<3i<n 
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It follows from [7J Lemma 5] that 

T =Q, 11 = 0, T 2 =Qo(IB 2 ) and T 3 = Q 0(£5 3 ). (2.2) 

We proceed to bound T s for s > 4. Let Y^' denote the sum over all sequences (ji, ■ ■ ■ E 

h,—,je 

{1, . . . , n} e with i distinct entries. Applying the multinomial theorem, we have 
where 



.mi, . . . , m£ 

miH \-me=s N 



k 



Let Ali be the set of all compositions m = (m 1; . . . , mg) of s such that = 1 for some i, 
and let M.2 be the set of all other compositions of s. For all m we have 

\B(m)\ < [n) e B s , 

using the falling factorial. For m E Mi, suppose as a representative case that mi — 1. Then 



3"f 

h,—>k 



= E' ^'•••Ct' E ft 
= - E' ^■■•C; 1 E 

h>—Jt-i ii£{h,—,h-\} 
where the last step uses the assumption 5^/=i /3j = 0. This shows that for m G .Mi we have 

\B(m)\ < £[n] e -iB s = 0(£/n)[n] e B s . 



Consequently 



s 

mi, ...,mt 



Furthermore 

s 



TOi, . . . , 777/ 



E 

Next, notice that for any fixed integer s > 4 

c -= E 



< 



.777-1 , . . . , 777/ 
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is the coefficient of x s in the Maclaurin expansion of s\ (e x — x) e . Since that expansion has 
nonnegative coefficients, C s < s\ri~ s (e v — rjY for any rj > 0. Substituting r\ = \J sji and 
using the fact that {ey® — y/x) 1 ^ < 2 for x > gives 

C s < s\ (s/£)- s/2 (eV^ - y/s/iy < s\ (2^/Ifs) s . 

Hence we have, for any fixed integer s > 4, 




Using ((H2D for s < 3 and (J2~3D for s > 4 > S ives 

Um = (fj (l + 0(n-^ 2 ) + 0(£/n) ^ ±B s t + 0(1) ^ B s (2^Ifs) s ) . 

Since B£ = o(l) and By/i = 0(n~ 1//4 ), the first sum in the above expression is 0(£/n) = 
0(n -1 / 2 ), while the second sum is at most 

J2B s £ s/2 = 0(B 4 £ 2 ) = din' 1 ). 

Hence 

U t (/3)=fy{l + d(n-W)), 

which matches the lemma for this range of I values. 

For the remaining range n—n 1 ^ 2 '^ < £ < n, we can apply the identity Ui({3) = U n -e(—(3), 
which is a consequence of V- f3j =0. The lemma is thus proved. □ 

2.2 Proof of the dense theorem (Theorem 11.41) 

Suppose that a, b > are constants such that a + b < |, and d is such that (11.21) holds 
and dj — d is uniformly 0(n 1//2+e ) for j = 1, . . . , n and some e > 0. In the following, we 
will assume that e is sufficiently small. Later in the proof we will infer that we can take 
e = e(a, b) for some function e(a, b) > 0, as required by Theorem 11.41 

Every vector z = (zi, . . . ,z n ) G {0, l} n is a potential diagonal of one of our matrices. 
Define \z\ — YTj=i z j anc ^ f° r ^ = 0, . . . , n let 

A t = {z e A : \z\= £}. (2.4) 
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If D£ and S have the same parity then 

G D (d,£) =^G(d- Dz). ( 2 - 5 ) 

We proceed by applying Theorem 11.31 to estimate G(d — Dz) and then summing the result 
over all z 6 Ag. Note that the average entry of d — Dz is d — D£/n. 

Let a be any constant such that a < a < | — b and let e = e (a,b) be the positive 
constant guaranteed by Theorem 11.31 Then for £ — 0, . . . , n we have 

. , , D£ , D£\ ^ n 
mm < a , n — a H > > 



n n J 3a log n 

for sufficiently large n. Provided e < Eq, we have that (dj — Dzj) — (d — D£/n) is uniformly 
O(n l ' 2+£0 ) for j = 1, . . . , n. So Theorem II .31 with the constants (a, 6) applies to every vector 
d — Dz, using the value Eq = £o(& ; b) guaranteed by that theorem. 

Next we will compare factors from the expression for G(d — Dz) given by (11. II) with 
corresponding factors from the formula for Gu(d,£) given in Theorem 11.41 Let denote 
the density of d — Dz for any z G Ag. That is, 

d D£ D(£-fx D n) 

n — 1 n[n — 1) n[n — 1) 

Also let 6j — dj — d for j — 1, . . . , n, which allows us to write R = YTj=i Then 



^™ /2 (l " ^)- D( - £)/2 exp ^ K ~^ n 2 + 0(n~ l 



= /in ' 1 - Ud v ;/ exp 

Using the expansion [m]^ = m fc exp ^— fc ^ m 1 ^ + 0(A; 3 /m 2 )j, valid when m — >■ oo such that 
k = o(m 2 / 3 ), we find that 

n-1 \ 

d 3- Dz i) _ „Dz jn „_\D(1- Zj )^f D(D - l)(ji D - Zj 



^ J (l-M^ j exp 



+ 0(n 



dj J rV.. N X D(,, -y.\2j\2 



D{fx D -z 3 )5 3 D^ D - Z j) 2 b 2 , ~^_ 3/2N 



\i D {\ - \i D )n 2/4,(1 - [i D ) 2 n 2 
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Since z 2 = Zj and YTj=i z j = ^ we obtain 

d i ) ( D(D-1) D(D - - 2fi D )(e - fji D n) DR 



2 2fi D (l -fi D )n 2(l-/i D ) 2 n 2 

^V/i D (l-AiD)n 2^(1-,, D )%V J 1 V 

Finally, apart from the 0(n~ b ) error term, the expression inside the exponential in (11. ip for 
d — Dz is 

1 (Ei=iK-^i-Mn-l)) 2 ) 2 1 ^ 



4 4A 2 (1-A,) 2 n 4 4 4/^(1 - /x D ) 2 ^ 

Combining these expressions gives 



+ 0{n- 1 ' 2 ). 



G{d -Dz) = A V{i) exp (3 jZj \ (2.6) 



where 



A = V2 (1 - to )'-"»)ft» +D " /2 f[ exp(0(,r») + O^)), (2.7) 

V(£) = fi" (1 - , D )<"-"^exp(i - - ^-l)(l-2^-^») 

\4 2 2/i£)(l — //d)tz 

L> 2 (£ - /i D n) 2 i? 2 



4/i D (l-/i D )n 2 2(l-/i D ) 2 n 2 4/^(1 - // D )V 
ft = _^ ^(1-2^ for j = l n 

Next we must sum over all z e 4^. Note that /3j = 0(n -1 / 2 ) for j = 1, . . . ,n, and the 
average of fii, . . . , fi n is 

= _ J*™* = S(n-). 
2/4(1 - /i D ) 2 ra 3 

Hence Lemma [2.11 applies and shows that 

5>p(f:/W) = (!) e X Ji ( 3 + i -^±^- ( 3r + d(n-^)) 
zeAt v i=i / V / V j=1 J 



n\ ( D 2 £(n-£)R D(l-2fi D )£R 



cxp 



2(j? D {l - Vo) 2 n 4 2^(l-/x D ) 2 n3 



+ 0(n- 1 / 2 )). (2.8) 



15 



Combining (EJ| and f|2TB|) — fT2T5]) gives 

G D (d,£) = A ^Vd /2 (1 - fi D ) {n - £)D/2 exp(g D (d,£) +0(n- 1 / 2 )) (2.9) 

where Q £>(d,£) is defined in the statement of Theorem II .41 for D G {1,2}. 

Next we will estimate Gd(d) by summing (12. 9p over allowable values of £. Recall the 
definition of £jj given in the theorem statement. Ignoring the factor A which is independent 
of £, we calculate 

£ [%T& - ^) D ^ /2 eMQn(d,£) + 0(n-^)) 

" / \ / \-DI/2 

= (1 -^) D " /2 £ J (tzM exp(g D (d,l D ) +0(^-4) +n-V2)) 
e=o \ / V 

n / \ / \ 

= (1 -M Dn/2 ex P (Q (d,4)) £ r) (^-j exp^Cn- 1 ^-^) +n- 1 / 2 )). 

(2.10) 

If \i — Id\ < n 1 ^ 2+v for some constant r] > then the error term in the corresponding 
summand is 0(n -1 / 2 ), so these summands are essentially terms from a binomial expansion. 
If \£-£ D \ > n 1/2+,? then 

" W f ' 1 ' » < exp(-fi(n 2r ')), exp^n" 1 ^-^))) = exp(0(l)), (2.11) 



.1 - V-D 

so the contribution from the tails of the sum is negligible. Therefore 

i i ex P - i 

I) \\ - [i D ) V V n 



£=0 



exp(5(„^)) YX ! X_^_) 



£=0 
D/2\n 



The preceding calculations hold for any sufficiently small e > 0, so in particular they 
hold for some e = e(a,b) such that e < Eq and the (^(n^ 1 ^ 2 ) error terms in (12. 7p . (12.91) and 
(I2.13P are all 0(n~ b ). Then the claimed formulae for Go(<i, £) follow immediately from (12 .7p 
and (El). 

Furthermore, multiplying ( 12.13ft by A (1 — jj,D) Dn exp(Qz)(<i, Id)) using ( 12.71) and sub- 
stituting D = 2 gives the desired formula for G2(d). 
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For D = 1, we must sum over only those values of I with the same parity as S. That 
is, we must replace (12. 12ft with a sum over just the even (or just the odd) values of I. By 
standard properties of the binomial distribution, the parity-restricted sum is half the full 
sum, within additive error 0(n~ b ), say. (This also follows from Lemma 13.31 when fii ^ |, 
and hence when jMi = 1/2 by analytic continuation.) The additive error can be absorbed 
into the relative error in (I2.13p . since the main factor there is Q(l). This gives the desired 
formula for Gi(d), completing the proof. □ 

3 The sparse case 

In this section we prove Theorem 11.51 
3.1 Some useful results 

First, we present two lemmas involving the Poisson binomial distribution, which we intro- 
duced in Section [L~T1 Let p = (pi, . . . ,p n ) satisfy < pi, . . . ,p n < 1 and let X be a random 
variable with distribution PB(p). The mean of X is X = E(X) = YTj=iPj- The following 
tail bounds are standard. 

Lemma 3.1. If X is a Poisson binomial random variable then, for any s >0, we have 

Pr(X - X < -s) < exp 

Pr(X - X > s) < exp^-X^ 
where <p(x) — (1 +x) log(l + x) — x. 

Proof. These bounds are attributed to Chernoff, see [TU| Theorems 2.1 and 2.8] and [SI 
Theorem 3.2]. □ 

Lemma 3.2. Let X be a random variable with Poisson binomial distribution PB(p) and 
mean X < n/(\ogn) 2 . For a fixed constant C > 0, let f : K. — >■ R be a function such that 

\f( x )\< C (-+W) for \x\<n. 

Then 

E(exp(/(X - X))) = 1 + E(f(X - X)) +0(E{f(X - X) 2 )) + rT n ^ n) 

= exp(E(f(X - X)) +0(E(f(X - X) 2 )) +n- n(logn) y (3.1) 
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In particular, the contribution to this expectation from values of X with \X — X\ > n 1 ^ 2 



is 



Pr(X = £)exp(f(£-X)) =n- Q{logn \ (3.2) 

Proof. Define g(x) = e^ x ' — 1 — f(x). Note that \g(x)\ < &^ x " for all x, which implies that 

\g(x)\ < exp( — h ) for \x\ < n. (3.3) 



n n 



We write 



where 



E(g(X-X)) = £ 1 + £ 2 + £ 3 



E 1= Yl P<X = l)g(l-X), 

i, \l-X\<nM 2 

£ 2 = P^(X = £)g(£-X), 

E 3 = J] Pr(X = £)g(£-X). 

1,1-XK-n 1 / 2 

In each of these sums, £ is a nonnegative integer in {0, . . . , n} which satisfies the additional 
constraint given. We now bound these three sums in turn. 

For Ei, note that when \£ - X\ < n 1 / 2 we have f{£ - X) < 2C = 0(1). Hence 

g{£-X)=0{f{£-X) 2 ) 

uniformly for all £ in this range. It follows that Yli w-xiKn 1 / 2 P r (^ = ■O /(^ ~~ -^) 2 — 
E(f(X -X) 2 ), and hence 

S 1 = 0(E(/(X-X) 2 )). (3.4) 

Now we consider £ 2 . Since X(p(s/X) is a decreasing function of X, and X < n/(\ogn) 2 
by assumption, Lemma 13.11 shows that 

Applying (I3.3P shows that £ 2 is bounded above by 

n max exp(Li(s)), 

n l/2 < s < n 

where 

T . . Cs 2 Cs n fs\og 2 n 

L u s ) = — + -in - r~^~ ^ 

n n 1 ' 2 log n V n, 



Now 

Tl ..2Cs C , / slog 2 n 
L[(s) = + - log 1 + 5_ 

which is negative for sufficiently large n for s G {n 1//2 ,n}. Also L'"(s) > for all s > 0, so 
it must be that L[(s) < for n 1//2 < s < n when n is sufficiently large. It follows that the 
maximum of L\ on the interval [n 1 ^ 2 , n] occurs at s — n l l 2 . Since Li(n 1 / 2 ) = — | log 2 n+0(l), 
we deduce that 

S 2 = n exp(-ft(log 2 n)) = rT n(logn) . (3.5) 



A bound on S3 can be obtained similarly. Using the first bound in Lemma I3.lt we find 

S 3 < n max exp(L 2 (s)), 

n 1 / 2 <s<n 

where 

T . . Cs 2 Cs s 2 log 2 n 
L 2 {s) = + -772 • 

By the same argument as before, the maximum of L 2 (s) occurs at s = n 1 / 2 for sufficiently 
large n, and we conclude that S3 = n~ n ( logn \ which together with (13.51) implies ( \3.2\j . 
Combining (13.21) and (13 .41) establishes (13.11) . completing the proof. □ 

For a given function / : {0, 1, . . . , n} — > K, define the polynomial / : M n — > M. by 

n 

f( Vl , ...,y n )= + •" + *») Hi/y (1 - Vj) 1 -**. 

(a;i,...,x n )e{0,l}" j=l 

(Note that this is indeed a polynomial in y\, . . . , y n , since 1 — Xj G {0, 1}.) In the case that 
< y%, ■ ■ ■ , y n < 1, we have 

f( yi ,...,y n )=E(f(Y)), (3.6) 

where Y is a random variable with distribution PB [(yi, ■ ■ ■ ,y n ))- 

The following lemma will be used when D — 1 to handle the parity restriction on the 
number of loops. 

Lemma 3.3. Fix (p\, . . . ,p n ) G [0, l] n snc/i that pj 7^ | /or j — 1, . . . ,n. Define 

rj = ~l-2 Pj 

for j = 1, ... ,n, and let 

n 

Z = JJ(l"2Pi)- 
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Then for p = 0, 1, 



E /(*! + ■■•+**) tlPs'^-Pi) 1 -' 

(x 1 ,...,x n )e{0,l} n j=l 
xiH \-x„=p (mod 2) 

= §/(pi,...,p„) + £ />i,..., r n ). 

Proof. Let X be a random variable with Poisson binomial distribution PB((pi, 
The probability generating function for X is 

P(w) = ^w'Pr(X = t) = -pj +pjw). 

t=0 j=l 

Note that 

n n 

f( Pl , ..., Pn ) = J2 f(t) Pr(X = t) = J2 /(*) [vf\P(w). 

t=0 t=0 

Now 

n n 

P(-w) = Y[(1-Pj- Pjw) = Z Yl (1 - Tj + r 3 w). 

3=1 3=1 

This expression has the same algebraic form as P(w), but with r 3 - in place of pj for 
1, . . . , n. Therefore, by comparison with (13.71) we have 



J2f(t) [w t ]P(-w) = Zf(r l ,...,r n ). 



t=o 

Hence we calculate that 



E f(xi+---+x n ) Hpfii-Pj) 1 ' 



(ae 1 ,...,!t n )6{0,l}» 3 = 1 

xi-\ \-x„=p (mod 2) 



E /(*) Pr ( x = *) 



t=0,...,n 
t=p (mod 2) 



E f(t)[w t ](iP(w) + (-iyiP(-w)) 



t — 0,...,n 

t=p (mod 2) 



E/w K](|PH + (-ir|p(-^)) 



t=0 



i/(p 1 ,...,p n ) + i(-l)^/(r 1 ,...,r n ), 



as claimed. 
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3.2 Proof of the sparse theorem (Theorem 11.51) 



We now prove Theorem 11.51 Assume throughout this section that 1 < d max = o(5' 1 / 3 ) and 
that S is even if D = 2. Furthermore, note that deleting vertices of degree zero does not 
affect either the value of Go(d) or the formulae for it given in Theorem 11.51 Hence we 
assume without loss of generality that dj > 1 for j = 1, . . . ,n. 
Let 



C C2 C2 C C4 C2 



2S AS 2 2S 4 AS 5 QS 3 



Using Stirling's approximation, Theorem 11.21 can be restated as follows: when S is even and 
1 < d max = o(5 1 /3), then 

G(d) = H(d)exp(0(d 3 m jS)), (3.8) 

uniformly as S — > oo. We proceed to estimate Go{d)/ H(d). 
Define 

A® = {ze {0, l} n : for j = 1, . . . , n, if dj < D then Zj = 0}, 
A® = {ze {0, l} n :\z\ = S (mod 2)} 



and let 

'A^nAW ifD = l, 
A® if D = 2. 

(Recall that \z\ denotes the number of entries of z equal to 1.) Then 



A 



G D (d) = H(d)J2 q ^ [ ^ 1 - (3-9) 



zeA 



Our strategy is to compare the ratio G(d — Dz)/H(d) to the ratio H(d — Dz)/H(d), which 
we now investigate. 



Lemma 3.4. For j = 1, . . . ,n, define 



a, 



^exp(4 + 7 ,), 



where 

D(D + 1) D(D + 2)S 2 DSl (£_ DS 2 , 
lj ~ 2S 2S 2 2S 3+ \S S 2 ' j 
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for j = 1, . . . , n, and A is defined by 

2^72 ~ Js? if D = 1, 



A 

Define K(z) by 



ifD = 2. 



H{d H(d) Z) =«■*(*(«)) IK- < 310 > 

Then there are functions K', K" : {0, 1, . . . , n} — > R which satisfy 

rp.jp. D 3f3 /J2 £2 M J3 \ 

Jftf), = -2H + ^ + ^ + O (^^f- + |j + (3.11) 

K'{\z\) < K{z) < K"{\z\) (3.12) 

for all z G A with \z\ < 5/3. 
Proof. Define the function 

M(d z) = _±Mf^L + — — ^ 2 ^^ 2 _|_ _ ^(z) 2 ^^) + Sf 5 3 



25i(z) 25 45i(z) 2 45 2 25i(2) 4 25 4 
2 (z) 4 SI S 3 (z) 2 S 2 



45i(z) 5 45 5 65i(z) 3 65 3 ' 
where S r (z) = Y^j=i[dj ~ Dzj] r for r = 1, 2, 3. Then 

H(d-Dz) ^/e\™ /2 / Dl \{S-DQ/2 n 



ff(d) exp(M( d , ,))(-) (l--) JJ <[*,]„ 

= exp(M( d ,z))e X p(^ + ^ + g))n( 



5-D/2 



Now 



5i(z) = S-Dl, 

S 2 (z) = S 2 - 2DW 1 + D(D + 1)£, 

S 3 (z) = S 3 - 3DW 2 + 3D(D + l)W l - D{D + \){D + 2)£, 
where W r = Y^j=i[dj]r %j for r = 1, 2. Making these substitutions gives 



'd 3 d 2 £ 2 
M (d,z)=0[^ + ^^- 
b b z 



) +J2^ z r 

' .7=1 



:r- 

Since the terms involving A cancel, this completes the proof. The lemma is in fact true for 
any A, but the value we have chosen will be useful in proving Lemma [3.61 □ 
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We now calculate some important quantities which will be needed later. 
Lemma 3.5. When D = 1, 

1 



.1 + 0* 



£>g(l + 

3=1 



s 



i 

2 



^2 2S*2 1S3 



C2 

^ 2 +0 



1 



5, 



So 5''? 

" +777-+V + 



245^2 25 2S 3 / 2 3S 3 / 2 2S 5 / 2 



C2 

^ 2 +0 



d 3 

"'max 

IT 

d 3 

ma: 

s 



When D = 2, 



St 



aj 



3=1 J 



S2 
S 



exp(0« ax /5)) ; 



^log(l + a J ) = ^exp(0(rfLx/^))- 
3=1 



Proof. For D = 1 we have 



Zi = 0(5- 1 / 2 ), 7, = o(%)+o(^)^, 



and find that 



n 1 



1 



+ 



2 85V2 2S 5 / 2 



C2 

^ 2 +0 



d 3 

max 

~s~ 



E2 °2 
fl i = 7F 



5*2 52 



7( 



j 5 5 3 / 2 



+ 1 



5-1/2 



5 



5'.') 3 So 1 
+ 7^+ + 



3=1 



53/2 53/2 5-1/2 

d 3 



5 



from which the result follows. When D = 2 we have 



j ~ s + v r ^ j ~ \ 52 ) 

3=1 3=1 



which imply the result in this case. 



(3.13) 



□ 



Next we calculate the sum of the right hand side of ( 13.101) over all z G {0, l} n (subject 



to a parity constraint if D = 1), after dividing by the factor nj=i(l + a j)- 
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Lemma 3.6. Let K* be either of the functions K', K" defined in Lemma 3.4 
If D = 1 then for p G {0,1}, 

E «p(*-(|,|)) ft t^- = * + 3^ + wr> + °(%)) 



^€{0,1}" i = i 

|z|=p (mod 2) 



IfD = 2 then 



n Zj 

ze{o,i} n j=i 3 



Proof. Define p = (pi, . . . , p n ) where = o 3 -/ (1+aj) for j = 1, . . . , n, and let X be a random 
variable with Poisson binomial distribution PB(p). Then 



n Zj 

a a 



exp(^(|z|)) Urfa- =H e MK*(X))). 
ze{o,i} n j=i i 

The expectation of X is X = Y2j=iPji which has been calculated for D = 1, 2 in Lemma [3751 
First suppose that D = 1. Recall from Lemma [3.51 that 

n 



From (13. 2 j) we know that 

E exp^dzl))!]^ 

l>3\/5 i =1 J 



n 2j 

' ; n -H(logn) 



|z|>3v / 5 

Next we observe that by Lemma 13.11 



n 2 

a,- 



For £ - y/S = O^S 1 / 3 ) we have, using lETTTj) . 

exp(JT(*)) = exp(-i + 0« ax /5))/(£), (3.14) 

where 

■' w on 09 o ' \ i e o ~ o ok/9 ~ 



32S 2 85 3 / 2 \16S 85 5 /2 485 3 / 2 

5 S 2 7 \ 41 5 2 5 



5 i/2 452 24S 32 8^ 3 / 2 485 1 /2 
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Since K*(£) = 0(1) for I = 0{VS), it follows that 

n Zj 



\z\=p (mod 2) 3=1 

-n(io g 



n 



")+exp(-I + 0« ax /5)) £ /(N)IIr+, 

\z\=p (mod 2) 3=1 J 



(3.15) 



We now apply Lemma 13.31 to estimate the sum on the right hand side. The small order 
moments of X are 



E(X 2 ) =X 2 + 



3=1 



E(X 3 ) 
E(X 4 ) 



X + 3X 2 + X 3 - 3 J2p 2 j ~ 3 ^ J2 p2 J + 2 E 



4 



3 = 1 



3=1 



3 = 1 



X + 7X 2 + 6X 3 + X 4 - (6X 2 - 18X + 7) J^P. 



3 = 1 



(n \ 2 n n 

Ep? + ( 8 *+ 12 )E^ 3 - 6 E 

3=1 J 3 = 1 3 = 1 



Substituting (13 . 1 3[) into these expressions gives 



E(X* 



JS-S 2 /S-\ + 0{dlJS^ 
S - 2S 2 /S^ + 0« ax ) 
S 3 / 2 - 3S 2 + p + 0(^^/2) 
^ S 2 + 4S 3 / 2 - 4S 1 / 2 S 2 + 0(dL*5') if = 4. 



if Jfc = 


0; 


if k = 


i; 


if Jfc = 


2; 


if fc = 


3; 


if fc = 


4. 



Hence 



/(Pi. • • • . P«) = E (/( X )) = ex P (^72- + + ( 



s 



(3.16) 



where / is the function obtained from / as in ( 13. 6ft . Let /& be the polynomial defined by 
= t k for all i G R, for = 1,2,3,4. Since E(X fc ) = /fc(pi, . . . ,p n ), replacing each pj 
with rj in (I3.16P leads to the following (abusing notation slightly to define the abbreviation 
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fi in the first line): 

n 

f 1 = f 1 (r 1 ,...,r n ) = Y,r j = 0(VS), 

n 

/ 2 (n, . . . ,r„) = A 2 + 5>(1 - r,) = 0(S), 



/sfa, ...,r„) = /i + 3/ 2 + ft - 3 J>f - 3/x E r i+ 2 E r ? = °^ /2 )' 

i=l j=i i=i 



• • • , r„) = /i + 7/7 + 6/f + ft - (6/7 - 18/! + 7) J] 



(n \ 2 7i 7i 

i=i 7 i=i i=i 



From this we conclude that /(n, . . . , r n ) = 0(1). Furthermore, 



Thus by Lemma [3.31 we obtain 

'dl 



e /(N>nT^M3^ + ^ +0 (- 

10 (mod 2) .7=1 /V V 



max 
IT 



|z|=p (mod 2) i= 

Combining this with (" 13 . 1 5[) establishes the lemma when D — 1. 

Next suppose that D = 2. Expanding Hf* around X gives 

K*(X) = h(X-X) + 0(d 3 m JS) (3.17) 
where ft, : R — >• R is a function which satisfies 

%)=o(%^y + o(iy (3.18) 

for \y\ < S. Recall our assumption that dj > 1 for j = 1, . . . , n, which implies that S > n. 
Hence the function h satisfies the conditions of Lemma [3.21 for some constant C > 0. We 
proceed to apply this lemma, specifically (13. ip . 

The second and fourth central moments of X are 



E((X-X) 2 ) =E>;(l-ft) = 0(X), 

n 

E((X - X) 4 ) = 3 E((X - X) 2 ) 2 + - Pi )(l - 6 Pj + 6p 2 ) = 0(X + X 2 ). 
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Recall from Lemma [3.51 that X = 0(d max ), and also note that \y\ <l + y 2 . From (I3.18P and 
flTTHD, we have 

E(h{X-X)) =0(rf max /5)(l+E((X-X) 2 )) =0{dlJS). 

Similarly, from (13 . 181) by applying (13. 191) and using the inequality {u + v ) 2 < 2{u 2 + v 2 ), we 
obtain 

E{h(X - X) 2 ) = 0(d 2 m JS 2 ) E((X - X) 2 ) + 0(S- 2 ) E((X - Xf) = 0« ax /S 2 ). 
Therefore ( 13. ip gives 

E(exp(h(X - X))) = exp(0(d 3 m JS)). 
This completes the proof when D = 2, using (I3.17p . □ 

We may now prove our main result in the sparse case. 

Proof of Theorem \1.5[ First suppose that S > n logn. Then Lemma E31 applies for all values 
of I. Furthermore, 

C = o(S - Dn) 

since S — Dn = Q{S), so (1 3 . 8 [) can be applied to d — Dz, for all z e A. Notice also that 
dj = whenever dj < D, so the sum of the right hand side of (13.101) over z £ A is equal to 
the sum over {0, 1}™ when D = 2, or over A^ when D — 1. Hence the result follows from 
(EES) using (|33J and Lemmas EaHSSl 

Now suppose that n < S < nlogn. We show that terms with \z\ > S/3 give a negligible 
contribution to Gz>(d). 

It is well known that when S is even, we can write 

rr, / n v -1 

G < d > = (s/2)M!H p(d) 

where -P(d) is a probability, and hence is at most 1. (Indeed, the exp(-) factor in (13. 8 p is 
an approximation to P(d) when d max = o(S 1 ^ 3 ), as proved in p2].) It follows by Stirling's 
approximation that 

S/2 / n 



G(d) = 0(l)(£) (j[d^ 



•3=1 

for any even value of S. Recall the definition of Ai from (12 .4p . For z e At we have 

/Q_ HA (S-D£)/2 ,n x-1 

G(d- J D 2 )=o(i)^^^j (n^-^) ! j • 
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Furthermore, 

5/2 n 



Hence 

^^) =0( l)exp( (0) (% £ ) D " 2 - 
Therefore, recalling that i < n and ignoring parity for an upper bound, 

E E^w^=o(d«p(o(0) E ©(W 2 

£=5/3 *S4« ^ ' (=S/3 ^ ' ^ ' 

= 0(l)exp(0(£J) E ("V"" 6 

£=5/3 ^ ' 

= 0(l)exp(0« ax ))2"S-^/ 18 

= 0(S^ (5) ). (3.20) 
Recall that (13. 8p applies when £ < S/3. Therefore, using (13. 8p and Lemma [3 A\ 
G D {d) 



iZYd) ^ Hid) 

5/3 „ 

= 0(S-^))+exp(OK a x/^)) £ E ex P(^W) IK 

£=0 Z6A« J=l 



Hence, by (j3J) , 



5/3 n 

0(^)) + exp(0(dL x /5)) E ex P(^'W) E II 

£=0 z£A e j=l 

~ H(d) 

5/3 n 

< 0(S- a &) +exp(0« ax /S)) E «p(**0) E IK ( 3 - 21 ) 

£=o ze4f i=i 

Next we would like to show that, in either the lower or upper bound in (13.21 j) . the sum over 
£ can be extended up to £ = n without affecting the answer significantly. Since every term 
is positive, zero is a lower bound for the tail of the sum. Again, we ignore the parity issue 
for an upper bound. Let K* be either K' or K". Firstly, note that since n < S we have 
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K*(£) = 0(£) uniformly for 5/3 < £ < n. Furthermore aj = o(S D//6 ) for j = l,...,n. 
Therefore 

n n n / \ 

E E II £ E ("j (e om s-""Y 

l=S/3 zeA e j=l e=s/3 ^ ' 



£=5/3 z£A e j=l l=S/3 

n 

<E C)s- D " 7 

£=S/3 



0(S-^ S) ) (3.22) 



as in (I3.20p . Combining this with (I3.2ip gives 



a J 



ze/1 j=i 

n 

< 0(S-^)+exp(0(£ v JS)) *M K "(\z\)) IR 

zeA j=i 

The result now follows from Lemma 13.51 and Lemma 13.61 □ 



4 Proof of Theorem 11.6 



Part (i). Under the conditions of Theorem 11.41 the distribution of Yd follows directly from 
(12.1 Op and (12. lip , noting in the case of D = 1 that the restriction of i to the same parity as 
S changes the normalizing factor by 2 to high precision, as explained in the last paragraph of 
Section I2T21 The formula for the expectation follows on summing £Prob( Yd = £), since the 
error term 0(e _ ™ n<1) ) contributes negligibly. To see that the same is true for the variance, it 
helps to use the cancellation-free formula 

Var(Z) = Prob(Z = k) Prob(Z = t) (k - I) 2 , (4.1) 
k<e 

which is true for all discrete random variables Z of finite variance (see for example [T2l p. 8]). 

Part (ii). Now suppose that the conditions of Theorem 11.51 hold and consider the case 
D — 1. Let X be a random variable with the Poisson binomial distribution PB(p), where 
Pj = ctj/ (1+dj) and a 3 - is defined in Lemma l3~4"l From (13.101) . (13.141) . we find that for I = \/S+ 
0(S 1 ^ 3 ), the distribution of Y\ is proportional to PB(p) to relative error 0(d 3 nax / S + S^ 1 ^ 3 ). 
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Moreover, the weight of both Y\ and X from \£ — y/~S\ > S 1 ^ 3 is e -5 " 11 ', and restriction 
of £ to the same parity as S contributes a factor of 2 to high precision as in the proof of 
Lemma 13.61 This gives, for £ — 0, . . . , n, 



Prob(F 1 = £) = (2 + 0(d 3 m JS + S^ 3 )) PB(p, t) + 0(e-^ (1) ). 

Next we show that the parameters p' in the theorem are sufficiently close to the parameters p. 
For each j, we find that 

p^expiOidlJS^ + S- 1 ))^. (4.2) 

By definition, 

pb(p,4= e (n^iK 1 -^))' 

where the sum is over subsets W C {1, 2, . . . , n} of size Applying (14. 2p . we find that 

PB(p,£) = exp(0(d 3 max /S + S^ 2 )) PB(p'J) 

for £ = O^). The tail past y/S is e" sn(1) for both PB(p) and PB(p'), by Lemma EH 
This completes the proof of the distribution for D — 1. 

The mean and variance follow as for part (i) to the same relative precision as the distri- 
bution, but we can do better by using the more accurate distribution analysed in the proof 
of Lemma 13.61 As we have shown in (13.141) . for £ = + 0(S 1 ^ 3 ), which excludes only 
exponentially small tails, 

Prob(F! = t) oc exp(0(d 3 ma JS)) Prob(A = £)f(£) (4.3) 

if £ has the same parity as 5*. Define the discrete random variable Z by 

Prob(Z = t) oc Prob(X = X + t)f(X + t), 

whenever X + t is an integer in [0, n] with the same parity as S; and Prob(Z = t) = 
otherwise. By (14. 3f) and the argument used in part (i) of this proof, 

E(Yi) = ex P (0(d 3 m JS))(X + E(Z)), 
Var(r 1 )=exp(0(dLx/5))Var(Z). 

For m > 0, define the central moment fi m = E((X — X)" 1 ) and the cumulant K m by 

oo 

log<f ) (t) = J2^t) m /m\, 

m=l 
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where (pit) = YYj=i{Pj elt + 1 — Pj) is the characteristic function of X. We find that n m = 
0(\/~S) for 2 < m < 6. Using the well-known expressions for the central moments in terms 
of the cumulants, and the explicit formulae (I3.19p . we find that 

^ = Vs- 2 -^-l + o(dLjs^), 
/i 3 = k 3 = o(Vs), 

/i 4 = 35 + 0(rf max 5 1 / 2 ), 

yU 5 = K 5 + 10k 3 K 2 = O(S), 

fi 6 = k 6 + 15k 4 k 2 + IOK3 + 15«ij = 0(S 3/2 ). 

Thus we calculate 

M = £'Prob(X = £) f(£) = ~ + * + J 2 ^ + 0« ax /S), 

n 

M x = ^'Prob(X = £) f(£) (£-X) = 0(dlJS + S^ 2 ), 

1=0 

M 2 = £ ' Prob(X = £) f(£) { £- X f = ^-^-\ + 0(dlJS^), 

1=0 

where the primes indicate that the sums are restricted to i having the same parity as S. The 
effect of the parity restriction is handled in the same way as in the proof of Lemma I3.6[ and 
in fact the first summation is equivalent to Lemma [3.61 Now we have that E(Z) = M\/Mq 
and Var(Z) = M 2 /M - E(Z) 2 . From these the mean and variance of Y\ follow. 

Finally we consider part (ii) in the case D = 2. Define X as before, with aj as in 
Lemma El By Lemma EH Prob(X > S 1 / 3 ) = 0{e~ snw ). The same bound holds 
for Pr(F 2 > S/3), using the argument leading to (13. 20 p . Combining this with (13. 2p and 
Lemma E31 shows that Pr(F 2 > 2S 1 ' 2 ) = 0(e~ sn(1) ). Finally, since K"(£) = 0(1) for 
£ = 0{S^ 2 ), we conclude that Pr(y 2 > S 1 ' 3 ) = 0{e~ sSlW ). 

Lemma [3.41 shows that for £ < S 1 ^ 3 , 

Prob(F 2 =i) = exp(0{d 3 ma JS + S" 1/3 )) Prob(X = £). 

By the argument above, the ratio of PB(p) to PB(p") for £ < S 1/3 is exp(0(d^ lflx / S 2/3 )) , 
since pj = exp(0(c/ 2 nax /5'))p" for all j. The given estimate of the distribution of Y 2 follows. 

To obtain the mean and variance of Y 2 , we use the sharper estimate 

Prob(F 2 = £) oc ex P (0(d 3 m JS)) Prob(X = £){! + £ 2 /S), 
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valid for t < dU&xS 1 / 4 by Lemma 13. 4[ with the weight of the tail £ > dU&xS 1 / 4 being expo- 
nentially small as usual. Using fl3TT6|) we find that E(X 2 ) = 0« ax ), E(X 3 ) = 0(d 2 nmx S 2 /S) 
and E(X 4 ) = 0(d 3 m&x S 2 /S), and so 

- ( d 2 \ 

VProb(X = £)(l + £ 2 /5) = l + o(=^), 

£ Prob(X = £)(! + £ = | + (%^) , 

£ Prob(X = £)(! + £ 2 /^) * = ff + | + O (%^) , 

and from these the expressions for E(la) and Var(Y"2) follow, recalling the cancellation-free 
variance formula 04. ip . □ 



5 A conjecture for regular graphs with loops 

In the case of D = 2 and d = (d,d, . . . ,d), an informal computation provides motivation for 
the sparse and dense enumeration formulae and suggests a more general conjecture. Since 
D = 2 we have d G {0, 1, . . . ,n + 1}. Recall the notations G 2 {n, d) = G2(d, d, . . . ,d) and 
fj, 2 — d/(n+ 1). 

Generate a random n-vertex graph by independently choosing each of the ( n ^ 1 ) possible 
edges (including loops) with probability /i 2 . Each d-regular graph has exactly nd/2 edges, 
so it occurs with probability 

^ /2 (l-V 2 f V )- nd/2 - (5-1) 
The event that a particular vertex has degree d has probability 

; ' lA t4 (i - *.r d +(T^ rfr 1 (i - ^r d+1 



d r Ay ^ ' \d-2 



n + l\ n- 1 d d+1 



din 



^{l-^) n - a+l . (5.2) 



If the vertex degrees were independent (which of course they are not), the number of graphs 
would be the n-th power of (15.21) divided by (15.11) . Noting that (1 — \/n) n — > e -1 , this gives 
a "naive" estimate 

G 2 (n,d) = e- 1 (*yX(tf(l-to) 1 -'»)W\ 
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We can see from Theorem 11.41 that Gzfo, d) is larger than G 2 (n, d) by a factor close to \/2e l l A 
whenever mm{d,n — d} > cn/ log n for some constant c > |. Less obviously, the same is 
true for 1 < d — o(n l l 2 ) by Theorem II .51 Recall that the same constant y/^e 1 ^ appears in a 
similar context for regular graphs without loops [16]. This leads us to investigate the region 
between the coverage of our sparse and dense theorems. 

Using the method described in [13], we computed the exact values of G 2 (n,d) for about 
150 nontrivial values of (n, d) up to n = 35. For example, 

G 2 (22, 10) = 7789744323722189254716829156528211234980743220762340514888. 

Numerical analysis of these values suggests the following analogue of [TBI Conj. 2]. 

Conjecture 1. Let d = d(n) satisfy 1 < d < n with nd even. Then 

uniformly as n — )■ oo ; where /i 2 — d/(n+ 1) and c = /i 2 (l — /i 2 )(n + 1). 

The numerical evidence suggests that in fact the term 0(n~ 2 ) always lies in the interval 
(-2/n 2 ,0) for n > 4. 
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